Effects of inhomogeneous activity of players and noise on cooperation in spatial public goods games 



O 
O 

(N 

> 
O 



Jian-Yue Guan^E Zhi-Xi Wu, 1 ' 2 '!! and Ying-Hai Wang 1 '! 

'institute of Theoretical Physics, Lanzhou University, Lanzhou Gansu, 730000, China 
2 Department of Electronic Engineering, City University of Hong Kong, Kowloon, Hong Kong, China 
(Dated: Received: date / Revised version: date) 

We study the public goods game in the noisy case by considering the players with inhomogeneous activity 
teaching on a square lattice. It is shown that the introduction of the inhomogeneous activity of teaching of the 
players can remarkably promote cooperation. By investigating the effects of noise on cooperative behavior in 
detail, we find that the variation of cooperator density pc with the noise parameter k displays several different 
behaviors: pc monotonically increases (decreases) with k; pc firstly increases (decreases) with k and then 
it decreases (increases) monotonically after reaching its maximum (minimum) value, which depends on the 
amount of the multiplication factor r, on whether the system is homogeneous or inhomogeneous, and on whether 
the adopted updating is synchronous or asynchronous. These results imply that the noise plays an important and 
nontrivial role in the evolution of cooperation. 
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I. INTRODUCTION 

Cooperation plays an important role in the real world, rang- 
ing from biological systems to economic and as well as so- 
cial systems IT]. Scientists from different fields of natural 
and social sciences often resort to evolutionary game theory 
llH as a common mathematical framework and the pris- 
oner's dilemma game (PDG) as a metaphor for studying co- 
operation between selfish and unrelated individuals yD. The 
PDG, which captures typical pairwise interactions in many 
cases, has become the leading paradigm to explain coopera- 
tive behaviorflllt]. As an alternative, the public goods game 
(PGG), which can be regarded as a PDG with more than two 
participators, attracts also much attention to study the emer- 
gence of cooperative behavior Jg]. In a typical PGG played 
by N agents J7Q, each of them must decide independently and 
simultaneously whether or not to invest to a common pool. 
The collected investment will be multiplied by a factor r, and 
then is redistributed uniformly among all players irrespective 
of their actual contributions. It was shown that for the values 
of 1 < r < N, the free riders or defectors (i.e., those who do 
not invest) will dominate the whole population lHIH]. 

To provide an escape hatch out of economic stalemate, 
Hauert et al. have introduced the voluntary participation in 
such public goods and found that it results in a substantial and 
persistent willingness to cooperate fich . Szabo et al. have 
studied the voluntary participation in PGG on square lattice 
10] and found that the introduction of loners leads to a cyclic 
dominance of the strategies and promotes substantial levels of 
cooperation, a remarkable increase of cooperation is also ob- 
served for those systems where the inhomogeneous imitation 
activity is introduced artificially to characterize the asymmet- 
ric and different influence of players to each other ifllll . In 
Ref. 11211 . Wu et al. have studied the PDG with the dynamic 



preferential rule and found that cooperation is substantially 
enhanced due to this simple selection mechanism. The effect 
of payoffs and noise on the maintenance of cooperative behav- 
ior has been studied in PDG where two types of players (A and 
B) with different activity of teaching are introduced on regular 
connectivity structures IU3I1 . It was found that the introduction 
of the inhomogeneous activity of teaching the players can re- 
markably promote cooperation Ultl . 

In the present, we study the effects of noise on the cooper- 
ative behavior by considering the players with the inhomoge- 
neous activity of teaching in PGG on the square lattice. Our 
goals are to investigate (i) whether the inhomogeneous activ- 
ity of teaching of the players can promote cooperation in PGG, 
and (ii) how the noise level affects the cooperative behavior 
in the PGG. Our observations suggest that, the cooperative be- 
havior is remarkably enhanced due to the introduction of the 
inhomogeneous activity of teaching. And we also find that the 
cooperator density pc changing with the noise level exhibits 
several different behaviors, which depends on the values of 
the multiplication factor r, on whether the system is homo- 
geneous or inhomogeneous, and also on whether the adopted 
updating is synchronous or asynchronous. 



II. THE MODEL AND SIMULATIONS 

Motivated by previous research work lfl2j[l3ll . we consider 
in our model two types of players (^4 and B) who are dis- 
tributed randomly on the sites of a square lattice before the 
start of the simulation. The fraction of players A and B are 
denoted by v and (1 — v), respectively. We consider the com- 
pulsory version of a spatial PGG (7], i.e., both types of players 
are only cooperators (C) or defectors (D). Each player inter- 
acts only with its four nearest neighbors (von Neumann neigh- 
borhood). Thus the size of interaction is g = 5. The payoff of 
each player depends on the number of cooperators n c in the 
neighborhood, i.e., the player x's payoff is 
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2 



Pa 



rn c 



0, if s x = D, 



(1) 



III. SIMULATION RESULTS AND ANALYSIS 



where s x denotes the strategy of the player x, and r represents 
the multiplication factor on the public good. The cooperative 
investments are normalized to unity and r > 1 must hold such 
that groups of cooperators are better off than groups of defec- 
tors. 

After each generation, the players try to maximize their in- 
dividual payoff by imitating (learning) one of the more suc- 
cessful neighboring strategies. Following previous studies 
JH 0, EH Ell H, player x will adopt the randomly chosen 
neighbor y's strategy s y with a probability depending on the 
payoff difference (p x — p y ) as 



1 



(2) 



where n denotes the amplitude of noise. The pre-f actor u> xy is 
given as 



= 1, 
to, 



if n y 
if n,, 
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where the value of u> (0 < uj < 1) characterizes the strength of 
reduced teaching activity if the site y is occupied by a player 
oftypeSlQ. 

Simulations were carried out for a population of N = 
100 x 100 individuals. We study the key quantity of cooperator 
density pc in the steady state. Initially, the two strategies of 
C and D are randomly distributed among the individuals with 
equal probability 1/2. The above model was simulated with 
both synchronous and asynchronous updating. Eventually, the 
system reaches a dynamic equilibrium state. The simulation 
results were obtained by averaging over the last 10 4 Monte 
Carlo (MC) time steps of the total 10 5 . Each data point results 
from an average over 20 realizations. 
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FTG. 1 : (Color online) The cooperator density pc vs the multiplica- 
tion factor r when v — 1 by using the synchronous and asynchronous 
updating. The noise parameter is k = 0.1. 



To facilitate comparison, we first consider the homoge- 
neous system (v = 1), i.e. the system for two types of players 
with the same teaching activity. Fig. Q] shows that the coop- 
erator density pc varies as the multiplication factor r when 
v = 1 by using synchronous and asynchronous updating. We 
can find from this figure that, below some threshold values of 
the multiplication factor (r < r s for synchronous updating 
and r < r as for asynchronous updating), the fraction of co- 
operators pc vanishes and the defectors reign; while for the 
values of multiplication factor r > r s and r > r as , the coop- 
erator density pc increases with r. For large r, pc can reach 
its maximum value 1, i.e., the system reaches the absorbing 
state of all cooperators, because large values of r are in favor 
of cooperators in the context of the PGG JUl- We compare 
the simulation results for synchronous and asynchronous up- 
dating and find that, when r is smaller than a certain value 
r c , pc is larger for asynchronous updating than that for syn- 
chronous updating, but the reverse is true for r > r c . Here 
we want to remark that in the process of our simulation, for 
synchronous updating, in each time step, all sites are updated 
simultaneously through competition with a randomly chosen 
neighbor. For asynchronous updating, however, in each time 
step, the individuals update their strategies one by one in a 
random sequence. It is known that in the context of the PGG, 
small values of r favor defectors and large r benefits cooper- 
ators Hlsl]. Comparing the average frequency of cooperators 
under the same values of r, we find that for small (large) r, 
defectors (cooperators) live better for synchronous updating 
than that for asynchronous updating. 
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FIG. 2: (Color online) The simulation results of cooperator density 
pc as a function of v for different values of uj when r = 4.8 and 
K = 0.1 by using synchronous (solid symbols) and asynchronous 
(hollow symbols) updating on the square lattice. The values of u> 
from top to bottom are respectively 0.01, 0.1, 0.2, and 1.0. 



In Fig. [2] we show the variation of cooperator density pc 
with the fraction v of individuals A for several different val- 
ues of u when r = 4.8 and k — 0.1 by using synchronous 
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(solid symbols) and asynchronous updating (hollow symbols). 
In our model, if u> = 1 the system is equivalent to a homo- 
geneous system (discussed in Fig. [TJ. For w < 1, i.e., the 
inhomogeneous system, the cooperator density pc increases 
monotonously until reaching the maximum value at v ~ 0.4 
(v ~ 0.33) for synchronous (asynchronous) updating, and 
then it decreases with v. And we also find that for the same 
value of v, the smaller the to is, the larger the pc will be. This 
phenomenon is similar to the results obtained in Ref. lfl3ll . 
We argue that it is a general phenomenon that the introduction 
of reducing activity of teaching can remarkably enhance the 
cooperative behavior in the context of whether PDG or PGG 
lfl5ll . From Fig. [2] we can also see that the qualitative results of 
the cooperative behavior are unchanged for both synchronous 
and asynchronous updating. Quantitatively, for the same value 
of v, the cooperator density pc is larger for synchronous up- 
dating than that for asynchronous updating when w < 1, but 
the reverse is true in the case of u> = 1. It indicates that, the 
enhanced level of cooperation due to the introduction of re- 
ducing activity of teaching depends on the detailed updating 
(synchronous or asynchronous). 
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FIG. 3: (Color online) The cooperator density pc vs the multipli- 
cation factor r in homogeneous (a),(b) and inhomogeneous (c),(d) 
systems by using the synchronous (a),(c) and asynchronous (b),(d) 
updating for three noise parameters: « = 0.01, 0.1, 1.0. 



In our simulations, we find that the level of noise, i.e., the 
magnitude of k, has a nontrivial effect on the evolution of co- 
operation. We show in Fig.[3jthe cooperator density pc vary- 
ing with the multiplication factor r in homogeneous [Figs. 
EJa) and (b)] and inhomogeneous [Figs. [2c) and (d)] sys- 
tems for several different values of k by using the synchronous 
[Figs. [3j a) and (c)] and asynchronous [Figs.[3jb) and (d)] up- 
dating. In the stationary state pc monotonously increases if r 
is increased for the same noise parameter. For the three noise 
parameter k = 0.01, 0.1, 1.0, in homogeneous systems [Fig. 
[3ja) and (b)], when r = 4.8, pc firstly increases and then de- 
creases with k, while for r = 5.1, pc monotonously decreases 



with k; for larger values of r = 6.0, the cooperator density 
displays different behaviors with k: pc monotonously de- 
creases with k for synchronous updating in Fig.[3ja); and pc 
firstly decreases and then increases with k for asynchronous 
updating in Fig. [3jb). In inhomogeneous systems (Fig. [3jc) 
and (d)), when r = 4.8, pc monotonously increases with k; 
while for r = 5.05, pc firstly decreases and then increases 
with k; for smaller values of r = 4.26, the cooperator density 
displays different behaviors with k: pc firstly increases and 
then decreases with k for synchronous updating in Fig.[3jc); 
and pc monotonously increases with k for asynchronous up- 
dating in Fig. [3jd). To investigate this in detail, we plot pc 
vs k for different values of r by using synchronous and asyn- 
chronous updating in homogeneous and inhomogeneous sys- 
tems respectively (see Fig.[U). 

Before making further discussions about our MC data, we 
want to give a brief description of the mean field result for the 
PGG. In a well mixed population (i.e., in the mean field case) 
with initial number of cooperator Npc, each cooperator and 
defector will get, respectively, rpc — 1 and rpc payoff. If a 
defector (cooperator) changes its strategy to be a cooperator 
(defector), then its payoff variation is r/N — 1 (and 1 — r/N). 
For any values of r larger than N, one can easily see that the 
transformation of defection to cooperation is preferred by the 
players, and the whole population will be dominated by coop- 
erators. In the reverse, for 1 < r < N, defectors dominate the 
population [5i, 8, 9]. Looking back on Fig. [3j it is clear that in 
the low noise limit k = 0.01, the MC data are in accordance 
with the results predicted by mean field theory JHH,[^], i.e., 
pc = if r < 5.0 and pc — 1 if r > 5.0 (N = g = 5 in the 
present case). 
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FIG. 4: (Color online) The cooperator density pc vs the noise param- 
eter n in homogeneous (a),(b) and inhomogeneous (c),(d) systems by 
using the synchronous (a),(c) and asynchronous (b),(d) updating for 
different values of r: the detailed values are shown in the plots. 



Figure [4] shows that the fraction of cooperators pc varies 
with the noise parameter k for different values of r in the 
two systems [the homogeneous in Figs. [2 a) and |2b), the 
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inhomogeneous in Figs. Uc) and Ud)] by using two updat- 
ing rules (the synchronous updating in Fig.|4](a) and (c), the 
asynchronous updating in Fig. H](b) and (d)). Let us first see 
the case in homogeneous systems. In Figs. [4] (a) and SJb), 
when r = 4.8, the cooperator density pc firstly increases 
monotonously until reaching the maximum value, and then 
decreases monotonously with k, which indicates that the co- 
operative behavior can be remarkably enhanced at the optimal 
level of noise both for synchronous and asynchronous updat- 
ing in homogeneous systems. For larger values of r, the con- 
centration of cooperators pc changing with k displays dif- 
ferent behaviors: It decreases monotonously with k for syn- 
chronous updating [see the curves for r = 5.1,5.5,6.0 in 
Fig. HI (a)]; while for asynchronous updating, when increas- 
ing the noise level, one can observe the valleys in the concen- 
tration of cooperators [see the curves for r = 5.1, 5.5, 6.0 in 
Fig. H](b)]. That is, a large level of noise evidently inhibits 
the persistence of cooperation for synchronous updating; but 
for asynchronous updating, the cooperative behavior can be 
inhibited mostly at certain value of k, which indicates that 
the cooperative behavior depends remarkably on the updating 
rule. In inhomogeneous systems, for synchronous updating 
[Fig. life)], the cooperator density pc varying with k displays 
three different behaviors: pc firstly increases with k and then 
decreases monotonously after reaching its maximum value for 
small value of r = 4.26; pc increases monotonously with 
K for moderate value of r = 4.8; pc firstly decreases with 
K and then increases monotonously after reaching its mini- 
mum value for large value of r = 5.05. This indicates that 
the changing behavior of pc with k are related closely to the 
multiplication factor r. But for asynchronous updating [Fig. 
EJd)], when r = 4.26, pc increases monotonously with k, 
which is different from the case for synchronous updating. 
Once again, this indicates definitely that the cooperative be- 
havior has strong dependence on the updating rule. 

In addition, for the synchronous updating and the same 
value of r, we can observe a peak in the concentration of 
cooperators when increasing the noise level in homogeneous 
system [see the curve for r = 4.8 in Fig.|4|a)], but the coop- 



erator density pc increases monotonously with the parameter 
K in inhomogeneous systems [see the curve for r = 4.8 in 
Fig. Etc)]. It indicates that the evolution of cooperation also 
depends on whether the system is homogeneous or inhomoge- 
neous. From the above behaviors, we can conclude that the ef- 
fects of noise on the cooperative behavior remarkably depend 
on the values of r, on whether the system is homogeneous or 
inhomogeneous, and on whether the updating is synchronous 
or asynchronous. 



IV. CONCLUSION 

In summary, we have studied the effects of noise on cooper- 
ation in spatial public goods games by considering two types 
of players with different activity of teaching. It was shown 
that the introduction of the inhomogeneous activity of teach- 
ing the players can remarkably enhance the persistence of co- 
operation. The introduction of the inhomogeneous activity of 
teaching can partially resolve the dilemma of cooperation and 
may shed new light on the evolution of cooperation in the so- 
ciety. By investigating the effects of noise on the coopera- 
tive behavior in detail, we found that the cooperator density 
pc varying with the noise level k displays several different 
behaviors, which remarkably depends on the values of r, on 
whether the system is homogeneous or inhomogeneous, and 
on whether the updating is synchronous or asynchronous. In- 
terestingly, when increasing the level of noise, we observed 
both peaks and valleys in the concentration of cooperators for 
the middle level of noise. Thus the effect of noise is a corre- 
lated factor in game dynamics, which plays an important role 
in the evolution of cooperation. 
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